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Several  algorithms  have  been  offered  to  track  the  Maximum  Power  Point 
when  we  have  one  maximum  power  point.  Moreover,  fuzzy  control  and 
neural  was  utilized  to  track  the  Maximum  Power  Point  when  we  have  multi - 
peaks  power  points.  In  this  paper,  we  will  propose  an  improved  Maximum 
Power  Point  tracking  method  for  the  photovoltaic  system  utilizing  a  modified 
PSO  algorithm.  The  main  advantage  of  the  method  is  the  decreasing  of  the 
steady  state  oscillation  (to  practically  zero)  once  the  Maximum  Power  Point 
is  located,  moreover,  the  proposed  method  has  the  ability  to  track  the 
Maximum  Power  Point  for  the  extreme  environmental  condition  that  cause 
the  presence  of  maximum  multi-power  points,  for  example,  partial  shading 
condition  and  large  fluctuations  of  insolation.  To  evaluate  the  effectiveness 
of  the  proposed  method,  MATLAB  simulations  are  carried  out  under  very 
challenging  circumstance,  namely  step  changes  in  irradiance,  step  changes  in 
load,  and  partial  shading  of  the  Photovoltaic  array.  Finally,  its  performance  is 
compared  with  the  perturbation  and  observation”  and  fuzzy  logic  results  for 
the  single  peak,  and  the  neural-fuzzy  control  results  for  the  multi-peaks. 
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1.  INTRODUCTION 

Solar  energy  is  alternative  source  of  energy  of  high  importance.  For  this  many  research  efforts  have 
been  concentrated  on  improving  the  efficiency  of  photovoltaic  (PV)  systems  and  accessibility  of  this 
technology.  The  PV  system  is  built  utilizing  one  or  more  panels,  the  power  converter  is  utilized  as  an 
interface  amongst  the  photovoltaic  array  and  the  load,  which  admits  a  Maximum  Power  Point  (MPP)  at  a 
precise  load  value.  This  last  is  not  stable  and  fixe,  it  subit  changes  as  parameters  change  the  comportment  of 
the  photovoltaic  system,  such  as  solar  radiation.  What  makes  it  essential  to  search  for  a  control  technique 
which  consists  in  adjusting  on  the  duty  cycle  of  the  converter  in  order  to  pursuit  the  MPP.  This  research 
presents  an  intelligent  approach  to  be  able  to  pursuit  the  Maximum  Power  Point  with  the  aid  of  technologies 
founded  on  particle  swarm  optimization  (PSO)  algorithm.  In  this  works  we  will  propose  an  improved 
maximum  power  point  tracking  (MPPT)  method  for  the  PV  system  using  a  modified  PSO  algorithm,  the 
results  got  from  this  algorithm  will  be  compared  with  the  Perturb  and  Observe  and  fuzzy  logic  algorithms  for 
the  single  peak  characteristics  of  the  P-V.  A  similar  comparison  is  done  for  the  multi-peak  characteristic  of 
the  P-V  between  the  PSO  and  neural-fuzzy  logic  algorithms  [1-2]. 
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2.  MODELING  OF  THE  PHOTOVOLTAIC  PANEL 
2.1.  Model  of  a  photovoltaic  cell 

The  equivalent  circuit  of  a  photovoltaic  cell  is  given  in  Figure  1  Maintaining  the  Integrity  of  the 
Specifications 
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Figure  1.  Equivalent  Circuit  of  a  Photovoltaic  Cell.  Figure  2.  An  Equivalent  Circuit  of  a  Photovoltaic  Field 


The  equivalent  circuit  of  Figure  1  is  composed  of  a  controlled  current  source  which  models  the 
effect  Photovoltaic  (the  current  generated  is  controlled  by  the  solar  illumination).  The  diode  represents  the 
effect  of  The  semiconductor  junction  of  the  cell  (models  use  two  diodes  in  parallel  for  more  precision).  Two 
resistors  (Rs,  Rp),  one  in  series  and  the  other  in  parallel,  produce  the  resistivity  effect  In  the  cell[l]. 

The  current  source  in  Figure  1  is  controlled  by  illumination  and  temperature: 


^pV  —  ^ 


(1) 


Or 

Ipv,n  is  the  nominal  current  (given  under  nominal  conditions:  T  =  25  °  C  and  G  =  1000W  /  m2),  KI 
is  the  Short-circuit  current  /  temperature  coefficient,  AT  =  T-Tn  (T  and  Tn  are  the  current  and  nominal 
temperature  Respectively),  G  and  Gn  are  the  current  and  nominal  illumination.  The  current  in  the  diode  is 
given  by: 
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With  10  the  saturation  current  of  the  diode: 
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Isc,n  is  the  nominal  short-circuit  current,  Voc,n  is  the  nominal  open-circuit  voltage.  KV  is  the  open 
circuit  voltage  /  temperature  coefficient,  a  is  a  constant  of  the  diode,  Vt  is  the  thermal  voltage  of  the  panel:  Vt 
=NskT/q,  where  Ns  is  the  number  of  cells  connected  in  series,  k  is  Boltzmann's  constant,  q  is  the  charge  of 
the  electron.  Rs  is  the  series  resistance  of  the  cell  (depends  on  the  material  used  to  construct  the  cell,  its  effect 
is  more  important  in  the  voltage  source  operating  region).  Rp  is  the  parallel  resistor,  its  effect  is  greater  in  the 
current  source  operating  region).  In  general,  the  values  of  the  resistors  are  not  given  by  the  constructor  and 
must  be  determined  by  iterative  algorithms  [3]. 

2.2.  Model  of  a  photovoltaic  filed 

A  photovoltaic  field  with  Npp  panel  in  parallel  and  Nss  panel  in  series  is  shown  in  Figure  2. 

The  output  current  is  then: 
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The  following  table  gives  the  parameters  of  the  photovoltaic  panel  used  in  the  following  paragraphs. 
Iscn  is  the  nominal  short-circuit  current  at  standard  test  conditions  (STC)  (25  °  C  and  1000  W  /  m2),  Ipvn  is 
the  nominal  current  of  the  panel  Ns  is  the  number  of  cells  (in  series)  for  each  panel,  a  is  a  constant  of  the 
diode,  Npp  is  the  number  of  modules  connected  in  parallel,  Nss  is  the  number  of  connected  modules  In 
series. 


Table  1.  PV  Panel  Settings. 


Parameter 

Value 

Iscn 

5.96 

Ipvn 

6.10 

Vocn 

64.2 

Ns 

96 

a 

1.3977 

Ki 

0.065%/°C 

Kv 

0.080V/°C 

Npp 

3 

Nss 

5 

Figure  3.  Configuring  the  PV  System. 


3.  Photovoltaic  Field  Control 

A  photovoltaic  field  is  the  combination  of  several  PV  panels  connected  in  series  and  parallel,  the 
slightest  defect  on  one  of  the  panels  can  cause  the  performance  degradation  of  the  entire  PV  field  [4-7]. 

That  is  why  research  has  proposed  an  architecture  or  a  controller  is  placed  at  each  set  of  small 
numbered  panels  in  order  to  manage  the  PV  field  decentrally.  Despite  the  effectiveness  of  the  proposed 
approach,  the  total  cost  of  such  an  installation  becomes  unaffordable.  Reason  one  considers  in  this  work  the 
architecture  or  a  single  control  must  manage  the  whole  field  PV. 

3.1  Configuring  the  Controller 

A  DC  /  DC  booster  converter  is  used  to  obtain  an  output  voltage  higher  than  that  provided  by  the  PV 
field.  The  PV  system  studied  is  illustrated  in  Figure  3.  Controller  parameters  are  given  in  Table  2. 


Table  2.  Controller  Settings. 


aComponent 

value 

Cin 

1  pF 

L 

0.005pH 

Cout 

24000  pF 

RL 

100  Q 
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Figure  4.  Characteristic  Curve  of  a  PV  Panel. 


3.2  Control  Algorithms 

1 .  P&O  algorithm: 

Several  algorithms  have  been  proposed  in  the  literature  to  extract  the  maximum  power  from  PV 
installations.  The  simplest  and  fastest  computational  time  algorithm  is  the  P  &  O  (Perturb  &  Observe) 
algorithm  which  is  based  on  the  disturbance  of  the  reference  current  to  visualize  the  direction  of  power 
progression.  Figure  4  illustrates  the  principle  of  this  control,  the  objective  is  to  reach  the  maximum  power 
point  (MPP)  in  a  minimum  of  time[l-2]. 

The  flowchart  of  the  'P  &  O'  command  is  given  in  Figure  5.  The  optimization  algorithm  requires  the 
presence  of  a  voltage  and  a  current  sensor  to  be  able  to  identify  the  direction  in  which  the  power  evolves. 


Figure  5.  Flowchart  of  the  P  &  O  optimization  algorithm. 


2.  Fuzzy  Logic 

Fuzzy  logic  asserts  itself  as  an  operational  technique.  Used  alongside  other  techniques  for  advanced 
control,  it  makes  a  discreet  entrance  but  popular  in  industrial  control  automation.  Fuzzy  controller  is  divided 
into  three  steps  illustrated  in  Figure  6. 


Figure  6.  Steps  of  fuzzy  control 


Figure  7.  Numeric  result  obtained  by  inference  block 
for  E=0  and  AE=0.5 
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Fuzzication:  This  section  transforms  inputs  into  fuzzy  variables.  The  inputs  are  transformed  into 
linguistic  variables  described  by  membership  functions.  These  functions  can  take  many  forms:  triangular, 
trapezoidal,  bell  or  other.  This  fuzzy  controller  uses  trapezoidal  and  triangular  membership  functions.  The 
input  and  output  variables  will  be  transformed  to  the  following  linguistic  variables:  1.  NG:  Grand  negative, 
2.  NM:  Middle  negative,  3.  N:  Negative,  4.  ZE:  Zero,  5.  P:  Positive,  6.  PM:  Middle  positive,  7.  PG:  Grand 
positive. 

Interference:  It  gives  the  relationship  between  input  and  output  variables  in  language  form.  There 
are  several  methods  of  inference:  MAX-MIN,  MAX-PROD  and  SUM-PROD.  Each  method  takes  its  own 
operators  and  each  rule  has  the  following  form: 

If  (Condition),  Then  (Conclusion) 

The  method  used  by  this  fuzzy  controller  is  the  MAX-MIN  method.  In  this  method  the  operators 
AND,  OR  and  Then  are  represented  by  the  MIN  and  MAX. 

Defuzzification:  During  this  step,  the  result  provided  by  the  fuzzy  inference  output  dD  is  converted 
to  real  value  or  numerical  value.  We  distinguish  several  methods  of  defuzzification:  the  maximum  method, 
the  weighted  average  method,  simplified  method  of  centroid,  the  centroid  method  and  others.  The  used 
method  by  the  used  fuzzy  controller  is  the  centroid  method.  This  method  will  then  convert  the  result  of  the 
fuzzy  output  variable  dD  to  the  duty  ratio  variations  that  will  be  acceptable  by  the  converter.  The  centroid 
method  calculates  the  center  of  gravity  of  the  fuzzy  surface  obtained  by  the  inference  block  (figure  7).  [9] 


4  APPLICATION  OF  PSO  ON  MPPT 

To  illustrate  the  application  of  the  PSO  algorithm  in  tracking  the  MPP  using  the  improved  technique 
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the  particle  position  xi  is  adjusted  where  the  velocity  component  vi  represents  the  step  size,  w  is  the  inertia 
weight,  Cl  and  C2  are  the  acceleration  coefficients,  rl,r2eU(0,l),  Pbesti  is  the  personal  best  position  of 
particle  ,  and  Gbesti  is  the  best  position  of  the  particles  in  the  entire  population.  If  position  is  defined  as  the 
actual  duty  cycle  while  velocity  shows  the  perturbation  in  the  present  duty  cycle,  then  can  be  rewritten: 

4+1 =<*>&  +c1rl\pLd  -4)+^(<4,d  -4) 

4=4+4+1  (6) 


In  the  case  of  PSO,  resulting  perturbation  in  the  present  duty  cycle  depends  on  Pbesti  and  Gbest.  If 
If  the  present  duty  cycle  (D)  is  far  from  these  two  duty  cycles,  the  resulting  change  in  the  duty  cycle  will  also 
be  large,  and  vice  versa.  In  the  last,  the  perturbation  in  the  D  is  always  fixed  but  in  PSO  it  varies  according  to 
the  position  of  the  particles.  With  proper  choice  of  control  parameters,  a  suitable  Maximum  power  point 
tracking  controller  using  PSO  can  be  easily  designed. 

4.1  MPPT  for  single  peak 

First  a  solution  vector  of  duty  cycles  with  Np  particles  is  determined,  i.e: 

dg=[dl,d2,d3,...,di] 

i=l,2,3,...,N  (7) 


where  dg  is  the  duty  cycle  of  the  DC -DC  convertor.  The  objective  function  is  defined  as: 

p(d^)>p(di]) 


(8) 


P(di  )  is  the  power  for  particle  i  at  iteration  k,  having  a  duty  cycle. 
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Figure  8.  PSO  particle  movements  in  searching  for 
the  MPP 


Figure  9.  MPPT  tracking  by  PSO  during  partial 
shading 


The  search  process  begins  to  transmits  three  duty  cycles  di  (i=  1,  2,  3)  to  the  power  converter.  In 
figure  8,  duty  cycles  dl,  d2,  and  d3  are  indicated,  with  triangular,  circular,  and  square  points,  respectively,  in 
the  first  iteration,  These  duty  cycles  served  as  the  Pbesti  .  Among  these,  d2  is  the  Gbest  that  offer  the  best 
fitness  value,  as  illustrated  by  Figure  8(a).  During  the  next  iteration,  the  resulting  velocity  is  only  due  to  the 
Gbest  term.  The  (Pbesti-d  (i))  factor  is  zero.  Moreover,  the  velocity  of  particle  Gbest  is  zero  due  to  the 
(Gbest-  d  (2))  factor  is  zero.  This  results  in  a  zero  velocity  and  accordingly  the  D  is  unchanged.  As  a  result, 
this  particle  will  not  effect  in  the  exploration  process.  To  avoid  such  situation,  a  small  perturbation  in  D  is 
allowed,  as  shown  in  Figure  8(b),  to  ensure  the  change  in  fitness  value.  Figure  8(c)  present  the  particles 
movement  in  the  third  iteration.  Due  to  the  fact  that  all  the  D  in  the  previous  iteration  attain  a  better  fitness 
value,  the  velocity  direction  of  these  particles  remains  unchanged  and  subsequently  they  move  toward  Gbest 
along  the  same  direction.  In  this  iteration,  all  duty  cycles  (di,i=  1,  2,  3)  arrive  at  Maximum  Power  Point  with 
a  low  value  of  velocity.  In  the  subsequent  iteration,  due  to  very  low  velocity,  the  value  of  the  D  is 
approaching  a  constant,  wherefore,  the  operating  point  will  be  keeped  and  the  oscillation  around  the  MPP 
diminishes. 

4.2  MPPT  for  multi-peaks 

When  the  photovoltaic  array  is  operating  in  a  uniform  solar  insolation,  the  resulting  P-V 
characteristic  curve  of  the  array  exhibits  a  single  maximum  power  point.  However,  under  partial  shading,  the 
P-V  curves  are  characterized  by  multiple  peaks,  i.e.,  with  several  local  and  one  global  peak  as  depicted  in 
figure  5.  In  this  example,  the  I-V  curve  is  characterized  by  four  stairs,  while  the  P-V  curve  is  characterized 
by  four  peaks.  The  latter  are  labeled  as  PI,  P2,  P3,  and  P4.  It  can  be  observed  that  the  time  derivative  of 
power  dP/dV  is  zero  for  the  global  as  well  as  all  the  local  peaks.  Furthermore,  the  slope  at  its  right  and  left 
sides  has  the  same  signs.  Since  all  conventional  MPPT  methods  are  based  on  the  slope  and  sign  value  of 
dP/dV,  the  algorithm  could  not  distinguish  the  local  (PI,  P2,  and  P4)  and  global  peaks  (P3)  correctly.  It  is 
very  likely  that  the  MPPT  is  being  forced  to  trap  into  the  local  peak,  resulting  in  reduced  output  power  and 
thus  greatly  deteriorates  the  efficiency  of  the  PV  system[5-8]. 


Figure  10.  MPPT  tracking  by  PSO  during  partial  shading 
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On  the  other  hand,  since  the  PSO  method  works  on  the  basis  of  search  technique,  the  global  peak 
can  be  tracked  without  any  difficulty.  Figure  9  depicts  the  tracking  capability  of  PSO  during  partial  shading. 
Similar  to  the  previous  P-V  curve  (see  figure  8,  the  proposed  method  transmits  three  duty  cycles,  which 
serve  as  Pbest  particles.  It  can  be  seen  that  the  voltage  and  current  contributing  to  these  initial  duty  cycles 
(Pbesti)  are  away  from  the  global  peak  (P3).  But  in  the  later  phases  of  iterations,  it  successfully  finds  the 
global  peak,  P3. 


5  RESULTS  AND  ANALYSIS 

Figure  10  shows  the  SIMULINK  model  of  the  proposed  system  in  single  peak,  figure  11  shows  the 
SIMULINK  model  of  the  proposed  system  in  multi -peaks. 


Figure  11.  Diagram  of  complete  system  of  PSO  in 
single  peak 


Figure  12.  Diagram  of  complete  system  of  PSO  in 
multi-peaks 


5.1  Comparaison  MPPT  for  single  peak 

Figures  12,  13  and  14  illustrate  the  three  models  used  for  the  different  simulations.  Figure  12 
presents  a  load  connected  to  PVG  with  the  aid  of  boost  type  DC -DC  converter  but  with  a  command  MPPT  of 
P  &  O  type.  Figure  13  presents  the  load  connected  to  the  same  PVG  but  this  time  using  the  fuzzy  command. 
Moreover,  figure  10  shows  the  MATLAB/Simulink  simulation  model  of  the  PV  system  used  with  the 
command  of  PSO  type.  SP  (solar  panel)  is  represented  by  a  Simulink  block,  and  DC -DC  converter  is  already 
described  above  in  figure  3.  The  frequency  of  the  impulsions  generator  chosen  is  40  kHz.  The  MPPT 
command  block  presents  the  used  command,  the  block  to  block  RC  MPPT  control  the  variation  of  the  duty 
cycle,  calculates  the  report  and  sends  it  to  HRV  function  of  “S -Function”  type,  which  admits  role  of  varying 
this  ratio  in  the  pulse  generator. 
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Figure  13.  Connection  with  the  P  &  O  algorithm 


Figure  14.  Connection  using  fuzzy  control 
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Figure  15.  The  maximum  tracked  power  of  the  PV 
model  using  PSO  algorithm 


Figure  16.  PSO,  P&O  and  fuzzy  logic  results  in 
tracking  the  MPP 


Figure  14  shows  the  variation  of  power  of  the  SP  as  a  function  of  time.  We  observe  the  effect  of  the 
applied  PSO  algorithm  on  the  maximum  power  point,  once  the  maximum  power  is  reached  the  system 
stabilized  at  this  value  with  an  oscillation  approximately  equal  to  8W.  We  will  compare  the  obtained  results 
of  the  PSO  algorithm  with  the  algorithm  "P&O"  and  the  algorithm  based  on  fuzzy  logic. 

Figure  15  shows  a  comparison  between  the  obtained  power  of  the  SP  of  PSO,  P&O  and  fuzzy  logic 
algorithms. 

From  the  results  obtained  in  figure  15  we  can  deduce  that  once  the  maximum  power  point  is 
reached,  the  systems  starts  oscillating  with  a  value  of  200  W  using  P&O  algorithm  and  50W  using  fuzzy 
logic  controller;  however  this  oscillation  decreases  to  8W  using  PSO.  We  can  deduce  that  in  tracking  the 
MPP  of  the  SP,  PSO  provides  high  accuracy  more  than  the  P&O  and  fuzzy  logic. 

5.2  MPPT  for  multi-peaks 

The  role  of  PSO  block  is  to  lead  the  system  to  the  MPPG  (maximum  power  point  global)  with 
margin  of  a  small  error,  and  then  stabilizes  the  system  at  the  maximum.  Figure  11  shows  the  diagram  of 
complete  system.  We  will  do  a  comparison  between  the  results  gone  using  the  PSO  algorithm,  and  the  results 
obtained  from  the  neural-fuzzy  logic  controller.  In  this  section,  we  present  two  different  simulations. 
Insolation  values  of  the  two  simulations  will  be  given  below. 

SIMULATION  1: 

Figure  4-13  shows  the  characteristic  curve  of  PV  system  under  the  following  insolation:  El  =  1000, 
E2  and  E3  =  700  =  300  W/m2 


Time(s) 


Figure  17.  P-V  characteristics  of  the  first  simulation 


Figure  18.  VPV  of  the  first  simulation  compared  to 
that  obtained  from  Neural  and  fuzzy  logic  control 
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Figure  19.  PPV  of  the  first  simulation  compared  to  Figure  20.  PPV(average)  of  the  first  simulation 

that  obtained  from  Neural  and  fuzzy  logic  control  compared  to  that  obtained  from  Neural  and  fuzzy 

logic  control 


SIMULATION  2: 

Figure  4-17  shows  the  characteristic  curve  of  photovoltaic  PV  system  under  the  following 
insolation:  El  =  750,  E2  =  950  and  E3  =  525  W/m2 


Figure  21.  P-V  characteristics  of  the  second  Figure  22.  VPV  of  the  second  simulation  compared  to 
simulation  that  obtained  from  Neural  and  fuzzy  logic  control 


Figure  23.  (a)  PPV  of  the  second  simulation  compared  to  that  obtained  from  Neural  and  fuzzy  logic  control, 
(b)  PPV(average)  of  the  second  simulation  compared  to  that  obtained  from  Neural  and  fuzzy  logic  control 


Due  to  non-uniform  solar  insolation,  the  resulting  P-V  characteristic  curve  of  the  SP  exhibits  multi¬ 
peaks  as  shown  in  figures  16  and  20.  Figures  17,  18,  21,  and  22  show  the  voltage  and  the  power  respectively 
of  the  PV  panel  obtained  using  the  PSO  algorithm  compared  to  that  obtained  from  the  neural-fuzzy  logic 
controller.  We  can  deduce  from  these  results  that  VPV  variation  using  the  neural -fuzzy  logic  controller  is 
about  40V  in  both  simulations,  where  it  is  decreased  by  a  considerable  value  using  the  PSO  algorithm  to  10V 
in  the  first  simulation  and  12V  in  the  second.  Moreover,  we  can  make  the  same  deduction  for  the  PPV  .Using 
neural-fuzzy  logic  controller,  the  power  oscillation  from  the  tracked  maximum  power  point  is  8000W  in  the 
first  simulation  and  400  W  in  the  second,  where  using  the  PSO  algorithm  it  is  300W  in  the  first  simulation 
and  200W  in  the  second.  Figures  20  and  23  show  the  variation  of  the  PPV(average)  using  PSO  algorithm  and 
neural-fuzzy  controller.  We  can  figure  out  that  once  the  maximum  power  point  (3000  W  in  first  simulation, 
4900  W  in  the  second  simulation)  is  reached  by  the  PSO,  it  is  stabilized  at  this  value;  however,  we  can  see  in 
neural-fuzzy  logic  controller,  once  the  maximum  power  point  is  reached,  the  system4 s  power  falls  to  a  value 
100W  less  than  the  maximum  power. 
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6  CONCLUSION 

This  paper  contains  an  approach  of  maximum  power  point  tracking  of  photovoltaic  systems  based 
on  particle  swarm  optimization  algorithm.  A  particle  swarm  optimization  (PSO)  with  the  capability  of  direct 
duty  cycle  is  used  to  track  the  MPP  of  PV  system.  A  PV  system,  designed  in  MATLAB  Simulink  was  used 
for  the  validation,  in  which  a  dc-dc  boost  converter  was  used  as  a  MPP  tracker.  The  performance  of  proposed 
controller  is  compared  with  “P&O”  and  fuzzy-logic  control  algorithms.  The  total  maximum  power  tracking 
using  the  PSO  algorithm  in  this  report  gives  a  better  result  than  the  “P&O”  and  fuzzy  ones  for  single -peak 
MPPT,  which  is  proved  by  simulation  experiments  in  this  report. 

In  the  case  of  partial  shading  of  PV  array  the  output  characteristics  presents  a  number  of  extreme 
points.  With  the  consideration  of  the  good  performance  of  the  PSO  in  multi -peak  function  optimization, 
global  optimization,  as  well  as  multi -variable  optimization,  the  PSO  algorithm  is  applied  to  the  MPPT  control 
algorithm.  The  used  PSO  algorithm  can  fast  and  accurately  realize  the  MPPT  in  PV  systems,  and  this  was 
proved  by  the  simulation  results  on  MATLAB  compared  to  that  obtained  by  the  neural -fuzzy  control 
algorithm.  However,  there  are  still  many  areas  for  improvement  for  the  PSO  algorithm  to  make  it  easy 
applicable  in  many  domains  especially  control. 
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